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Abstract. The toric Hilbert scheme parametrizes all algebras isomor- 
phic to a given semigroup algebra as a multigraded vectorspace. All 
components of the scheme are toric varieties, and among them, there is 
a fairly well understood coherent component. However, it is unknown 
whether toric Hilbert schemes are always connected. In this chapter 
we illustrate the use of MacauJay 2 for exploring the structure of toric 
Hilbert schemes. In the process we will encounter algorithms from com- 
mutative algebra, algebraic geometry, polyhedral theory and geometric 
combinatorics. 



Introduction 

Consider the multigrading of the polynomial ring R = C[x\, . . . , x n ] spec- 
ified by a non-negative integer d x n-matrix A = (ai, . . . ,a n ) such that 
degree(xi) = Oj E N d . This defines a decomposition R = ®beNA^b, where 
NA is the subsemigroup of N d spanned by ax, ■ ■ ■ ,a n , and Rf, is the C-span 

of all monomials x u = x™ 1 ■ ■ ■ x^ 1 with degree Au = axux H + a n u n = b. 

The toric Hilbert scheme HilbA parametrizes all A-homogeneous ideals 
I C R with the property that (R/I)b is a 1-dimensional C-vectorspace, 
for all b G NA We call such an ideal / an A- graded ideal. Equivalently, I 
is ^4-graded if it is ^4-homogeneous and R/I is isomorphic as a multigraded 
vectorspace to the semigroup algebra C[Nt4] = R/I a, where 

I A ■= ( x u -x» : Au = Av) C R 

is the toric ideal of A. It follows from the definition that ^4-graded ideals 
are generated by monomials and binomials in R. 

Our running example throughout this chapter is the following 2 x 5-matrix: 

m A f 1 1 1 1 A 

[L) A ~ [p 1 2 7 8J 

which we input to Macaulay 2 as a list of lists of integers. 

11 : A = {{1,1, 1,1,1}, {0,1,2,7,8}}; 

The toric ideal of A lives in the multigraded ring R = C[a, b, c, d, e]. 

12 : R = QQ [a. . e ,Degrees=>transpose A]; 

13 : describe R 

o3 = QQ [a, b, c, d, e, Degrees => {{1, 0}, {1, 1}, {1, 2}, {1, 7}, {1 • • • 
o3 : Adjacent 

1 
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We use Algorithm 12.3 in p(J to compute Ia- The first step is to find a 
matrix B whose rows generate the lattice ker%(A) := {i £ Z* : Ax = 0}. 

i4 : B = transpose syz matrix A 

o4 = I 1 -2 1 1 
10 5 -6101 
10 6 -7 11 

3 5 
o4 : Matrix ZZ < ZZ 

Although in theory any basis of kerz(A) will suffice, in practice it is more 



efficient to use a reduced basis [18, §6.2] which can be computed using the 
basis reduction package LLL.m2, a built-in feature of Macaulay 2 . The ad- 
vantage of a reduced basis may not be apparent in small examples. However, 
as the size of A increases, it becomes increasingly important for the termi- 
nation of Algorithm 12.3. (To appreciate this, consider the matrix (R) from 
Section 4.) 

15 : load "LLL.m2"; 

16 : B = transpose matrix LLL syz matrix A 

06 = I 1 -1-11 I 
11-10 -11 I 
12 -3 2 -1| 

3 5 
06 : Matrix ZZ < ZZ 

A row b = b + — b~ of B is then coded as the binomial x b+ — x b G R, and 
we let J be the ideal generated by all such binomials. 

±7 : toBinomial = (b,R) -> ( 

top := 1_R; bottom := 1_R; 

scan(#b, i -> if b_i > then top = top * R_i~(b_i) 

else if b_i < then bottom = bottom * R_i"(-b_i)); 
top - bottom) ; 

18 : J = ideal apply (entries B, b -> toBinomial (b ,R) ) 

2 2 3 

08 = ideal (- c*d + b*e, - b*d + a*e, ad - c e) 
08 : Ideal of R 

The toric ideal equals (J : (xi • • • x n )°°) which is computed via n successive 
saturations as follows. 

19 : scan(gens ring J, f -> (J = saturate (J,f) ;) ) 

Putting the above pieces of code together, we get the following procedure 
for computing the toric ideal of a matrix A. 

110 : toricldeal = (A) -> ( 

n := #(A_0); 

R = QQ [vars(0 . .n-1) ,Degrees=>transpose A,MonomialSize=>16] ; 
B := transpose matrix LLL syz matrix A; 
J := ideal apply (entries B, b -> toBinomial (b ,R) ) ; 
scan(gens ring J, f -> (J = saturate( J , f ) ; ) ) ; 

); 

In our example, Ia = {cd — be, bd — ae, b 2 — ac, a 2 d 2 — c 3 e, c 4 — a 3 e, be 3 — 
a 3 d,ad 4 — c 2 e 3 ,d 6 — ce 5 ) which we now compute using this procedure. 

111 : I = toricldeal A 

2 2 2 3 4 3 

oil = ideal (c*d - b*e, b*d - a*e, b - a*c, ad - ce, c - ae, b*c • • • 

oil : Ideal of R 
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This ideal defines an embedding of P as a degree 8 curve into P 4 . We 
will see in Section 3 that its toric Hilbert scheme HilbA has a non-reduced 
component. 

We recommend |2(| as an introductory reference for the topics in this 
chapter. The study of toric Hilbert schemes for d = 1 goes back to Arnold 
[|lj and Korkina et.al.[[ll]], and it was further developed by Sturmfels (|i~9|] 
and [2C, §10]). Peeva and Stillman [15| introduced the scheme structure 
which gives the toric Hilbert scheme its universal property, and from this 
they derive a formula for the tangent space of a point on HilbA- Maclagan 
recently showed that the quadratic binomials in M, §5] define the same 
scheme as the determinantal equations in [15]. Both of these systems of 
global equations are generally much too big for practical computations. In- 
stead, most of our algorithms are based on the local equations given by 
Peeva and Stillman in [14] and the combinatorial approach of Maclagan and 
Thomas in gg|. 

This chapter is organized into four sections and two appendices as fol- 
lows. The main goal in Section 1 is to describe an algorithm for generating 
all monomial A-graded ideals for a given A. These monomial ideals are 
the vertices of the flip graph of A whose connectivity is equivalent to the 
connectivity of HilbA- We describe how all neighbors of a given vertex of 
this graph can be calculated. In Section 2, we explain the role of polyhedral 
geometry in the study of HilbA- Our first algorithm tests for coherence in 
a monomial A-graded ideal. We then show how to compute the polyhedral 
complexes supporting A-graded ideals which in turn relates the flip graph 
of A to the Baues graph of A. For unimodular matrices, these two graphs 
coincide and hence our method of computing the flip graph can be used to 
compute the Baues graph. Section 3 explores the components of HilbA via 
local equations around the torus fixed points of the scheme. We include 
a combinatorial interpretation of these local equations from the point of 
view of integer programming. The scheme HilbA has a coherent component 
which is examined in detail in Section 4. We prove that this component 
is, in general, not normal and that its normalization is the toric variety of 
the Grobner fan of I a- We conclude the chapter with two appendices, each 
containing one large piece of Macaulay 2 code that we use in this chapter. 
In Appendix 1 we describe polarCone .m2 which is a procedure to convert a 
generator representation of a polyhedron to an inequality representation and 
vice versa. In Appendix 2 we explain a simplified version of the procedure 
minPres .m2, available in Macaulay 2 , for computing minimal presentations 
of polynomial quotient rings. The main ingredient of this package is the 
subroutine removeRedundantVariables which is what we use in this chap- 
ter. 
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1. Generating monomial ideals 

We start out by computing the Graver basis GrA, which is the set of 
binomials in Ia which are minimal with respect to the partial order 

x u — x v < x u — x v : <^=^ x u divides x u ' and x v divides x v ' . 

The set Gr& is a universal Grobner basis of Ia and has its origins in the 
theory of integer programming ||. It can be computed using [2C, Algorithm 
7.2], a Macaulay 2 version of which is given below. 

112 : graver = (I) -> ( 

R := ring I; 

kk := coef f icientRing R; 
n : = numgens R ; 

— construct new ring S with 2n variables 
S := kk[Variables=>2*n,MonomialSize=>16] ; 
toS := map(S,R, (vars S)_{0. .n-1}) ; 

toR := map(R,S,vars R I matrix(R, {toList (n: 1)}) ) ; 

— embed I in S 
J := gens toS I; 

— construct the toric ideal of the Lawrence 

— lifting of A 
i := 0; 

while i < n do ( 

wts := join(toList(i:0) ,{1} ,toList (n-i-1 : 0) ) ; 

wts = join(wts.wts) ; 

J = homogenize(J,S_(n+i) ,wts) ; 

i=i+l; 

); 

J = ideal J; 

scan(gens ring J, f -> (J = saturate (J,f) ;)) ; 

— apply the map toR to the minimal generators of J 

J = matrix entries toR mingens J; 

p := sortColumns J; 

J_p) ; 

The above piece of code first constructs a new polynomial ring S in n more 
variables than R. Assume S = C[x\, . . . ,x n ,yi,... ,y n ]. The inclusion map 
toS : R —* S embeds the toric ideal I in S and collects its generators in 
the matrix J. A binomial x a — x b lies in GrA if and only if x a y b — x b y a 
is a minimal generator of the toric ideal in S of the (d + n) x 2n matrix 

A{A) := ( T t ) which is called the Lawrence lifting of A. Since u £ 

kerz(A) 44> (u, — u) £ ker%(A(A)), we use the while loop to homogenize 
the binomials in J with respect to A(A), using the n new variables in S. 
This converts a binomial x a — x b G J to the binomial x a y b — x b y a . The ideal 
generated by these new binomials is also labeled J. As before, we can now 
successively saturate J to get the toric ideal of A(A) in S. The image of 
the minimal generators of this toric ideal under the map toR : S — > R such 
that and yi i— > 1 is precisely the Graver basis GrA- This list is then 

sorted and output. 

In our example GrA consists of 42 binomials. 

113 : Graver = graver I 

ol3 = {0, 0} I -cd+be -bd+ae -b2+ac -cd2+ae2 -a2d2+c3e -c4+a2bd -c4+a3 • • • 

1 42 
ol3 : Matrix R < R 

Returning to the general case, an element b of NA is called a Graver 
degree if there exists a binomial x u — x v in the Graver basis GrA such that 
Au = Av = b. If b is a Graver degree then the set of monomials in R\, 



ALGORITHMS FOR THE TORIC HILBERT SCHEME 



5 



is the corresponding Graver fiber. In our running example there are 37 
Graver fibers each corresponding to a different Graver degree. To compute 
the Graver fibers of A, we first construct the Doubleldeal DI := (x a x b : 
x a — x b 6 Gr a) which is contained in every monomial ideal of HzWa- Since 
our purpose in constructing Graver fibers is to choose standard monomials 
from each fiber, it suffices to list the monomials in each Graver degree that do 
not lie in DI. Since R is multigraded by A, we can obtain such a presentation 
of a Graver fiber by simply asking for the basis of R in degree b modulo DI. 

114 : graverFibers = (Graver) -> ( 

Doubleldeal := (I) -> ( trim ideal ( 
apply (numgens I, a -> ( 

f := I_a; leadTerm f * (f - leadTerm f))))); 
DI := Doubleldeal ideal Graver; 
R := ring Graver; 
result := new MutableHashTable ; 
scan (degrees source Graver, d -> ( 
if not result#?d then 

result#d = compress (basis(d,R) DI))); 

result) ; 

115 : fibers = graverFibers Graver; 

116 : peek fibers 



ol6 = MutableHashTable{{2 
{2 
{2 
{3 
{4 
{4 
{4 
{5 
{5 
{5 
{5 
{5 
{5 
{5 
{5 
{6 
{6 
{6 
{6 
{6 
{6 
{6 
{6 
{6 
{6 
{6 
{6 
{7 
{7 
17 
{7 
{7 
{7 
{7 
{7 
{8 
{8 



2} = 
8} = 
9} = 
16} 
14} 
7} = 
8} = 
10} 
14} 
16} 
21} 
22} 
28} 
7} = 
8} = 
12} 
14} 
18} 
21} 
24} 
28} 
30} 
35} 
36} 
42} 
7} = 
8} = 
14} 
21} 
28} 
35} 
42} 
49} 
7} = 
8} = 
56} 
8} = 



-> {0, 
-> {0, 
-> {0, 
=> {0, 
=> {0, 

■■> {0, 
■■> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0.. 
-> {0, 
-> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0.. 

=> {0, 

■■> {0, 

-> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

=> {0, 

■■> {0, 

■■> {0, 
=> {0, 
-> {0, 



0} I 
0} I 
0} I 
0} 
0} 
0} I 
0} I 
0} 
0} 
0} 
0} 
0} 
0} 
0} I 
0} I 
0} 
0} 
0} 
0} 
0} 
0} 
0} 
0} 
0} 
0} 
0} I 
0} I 
0} 
0} 
0} 
0} 
0} 
0} 
0} I 
0} I 
0} 
0} I 



ac b2 I 
ae bd I 
be cd I 

I ae2 bde cd2 I 
I a2d2 c3e I 

a3d bc3 I 

a3e a2bd c4 I 

I a3ce a2b2e a2bcd ab3d c5 I 
I a3d2 ac3e b2c2e bc3d I 
I a3e2 a2cd2 ab2d2 c4e I 
I a2d3 bc2e2 c3de I 
I a2d2e abd3 c3e2 I 
I ad4 c2e3 I 
a4d abc3 b3c2 I 
a4e a3bd ac4 b2c3 I 
a3c2e a2bc2d ab4e b5d c6 I 
a4d2 a2c3e abc3d b4ce b3c2d 
a3ce2 a2b2e2 a2c2d2 b4d2 c5 
a3d3 abc2e2 ac3de b3ce2 bc3 
a3e3 a2cd2e abcd3 b3d3 c4e2 
a2d4 ac2e3 b2ce3 c3d2e I 
a2d2e2 acd4 b2d4 c3e3 I 
ad5 bce4 c2de3 I 
ad4e bd5 c2e4 I 
ce5 d6 I 
a5d a2bc3 b5c I 
a5e a4bd a2c4 b4c2 I 
a5d2 a3c3e a2bc3d b6e b5cd 
a4d3 a2bc2e2 a2c3de abc3d2 
a3d4 a2c2e3 ac3d2e b4e3 bc3 
a2d5 abce4 ac2de3 b3e4 c3d3 
ace5 ad6 b2e5 c2d2e3 I 
be6 cde5 d7 I 
a6d a3bc3 b7 I 
a6e a5bd a3c4 b6c I 
I ae7 bde6 cd2e5 d8 I 
a7e a6bd a4c4 b8 I 



ol6 : Met 

For example, the Graver degree (8, 8) corresponds to the Graver fiber 



r 7 6ij 4 4 
| a e , a bd , a c 



a 3 6 2 c 3 , 



ab 6 c, &!}. 



Our Macaulay 2 code only outputs the four underlined monomials, in the 
format | a7e a6bd a4c4 b8 I . The three non-underlined monomials lie in 
the Doubleldeal. Graver degrees are important because of the following 
result. 
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Lemma 1.1. [20, Lemma 10.5]. The multidegree of any minimal generator 
of any ideal I in HUoa is a Graver degree. 

The next step in constructing the toric Hilbert scheme is to compute all its 
fixed points with respect to the scaling action of the n-dimensional algebraic 
torus (C*) n . These fixed points are the monomial ideals M lying on HilbA- 
Every term order -< on the polynomial ring R gives such a monomial ideal: 
M = m^(I^), the initial ideal of the toric ideal I a with respect to -<. We 
call these the coherent monomial ideals on HilbA- They can be computed 
by |2C, Algorithm 3.6] applied to I a- A refinement and fast implementation 



can be found in the software package TiGERS by Huber and Thomas p0[ . 

Now we wish to compute all monomial ideals M on HilbA regardless of 
whether M is coherent or not. For this we use the procedure generateAmonos 
given below. This procedure takes in the Graver basis GrA and records the 
numerator of the Hilbert series of I a m trueHS. It then computes the Graver 
fibers of A, sorts them and calls the subroutine selectStandard to generate 
a candidate for a monomial ideal on HilbA- 

117 : generateAmonos = (Graver) -> ( 

trueHS = poincare coker Graver; 

fibers = graverFibers Graver; 

fibers = apply(sort pairs fibers, last); 

monos = {}; 

selectStandard(f ibers , ideal (0_ (ring Graver))); 

The input to the subroutine selectStandard are the Graver fibers given 
as a list of matrices and a monomial ideal J which should be included in 
every ^-graded ideal that we generate. The subroutine then loops through 
each Graver fiber, and at each step selects a standard monomial from that 
fiber and updates the ideal J by adding the other monomials in this fiber to 
J. The final J output by the subroutine is the candidate ideal that is sent 
back to generateAmonos. It is stored by the program if its Hilbert series 
agrees with that of I a- 

118 : selectStandard = (fibers, J) -> ( 

if #f ibers == then ( 

if trueHS == poincare coker gens J 

then (monos = append (monos , flatten entries mingens J)); 
) else ( 

P := fibers_0; 
fibers = drop(f ibers , 1) ; 
P = compress(P '/, J); 
nP := numgens source P; 

— nP is the number of monomials not in J. 
if nP > then ( 

if nP == 1 then selectStandard(f ibers, J) 
else ( — remove one monomial from P.take the rest. 
P = flatten entries P; 
scan(#P, i -> ( 

Jl := J + ideal drop(P , {i , i}) ; 
selectStandard(f ibers , Jl))))); 

)); 

All the monomial A-graded ideals are stored in the list monos. Below, we 
ask Macaulay 2 for the cardinality of monos and its first ten elements. 

119 : generateAmonos (Graver) ; 

120 : #monos 
o20 = 281 
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i21 : scan(0..9, i -> print toString monos#i) 

{c*d, b*d, b~2, c~3*e, c~4, b*c~3, c~2*e~3, b*c~2*e~2, b*c*e~4, d~6} 

{c*d, b*d, b~2, c~3*e, c~4, b*c~3, c~2*e~3, b*c~2*e~2, c*e~5, b*c*e~4, • • • 

{c*d, b*d, b~2, c~3*e, c~4, b*c~3, c~2*e~3, b*c~2*e~2, c*e~5, b*c*e~4, • • • 

{c*d, b*d, b~2, c~3*e, c~4, b*c~3, c~2*e~3, b*c~2*e~2, c*e~5, b*c*e"4, • • • 

{c*d, b*d, b~2, c~3*e, c~4, b*c~3, c~2*e~3, b*c~2*e~2, d"6, a*d~5} 

{c*d, b*d, b~2, c~3*e, c~4, b*c"3, b*c~2*e~2, a*d"4, d"6} 

{c*d, b*d, b~2, c~3*e, c~4, b*c"3, a*d"4, a~2*d~3, d~6} 

{c*d, b*d, b"2, a"2*d"2, c"4, b*c"3, a*d"4, d"6} 

{c*d, b*d, b"2, a"2*d"2, a"3*d, c~4, a*d"4, d"6} 

{c*d, b*d, b"2, a"3*e, a"2*d"2, a"3*d, a*d"4, d~6} 

The monomial ideals (torus-fixed points) on HilbA form the vertices of the 
flip graph of A whose edges correspond to the torus-fixed curves on HilbA- 
This graph was introduced in [12] and provides structural information about 



HilbA- The edges emanating from a monomial ideal M can be constructed as 
follows. For any minimal generator x u of M, let x v be the unique monomial 
with x v ^ M and Au = Av. Form the wall ideal which is generated by 
x u — x v and all minimal generators of M other than x u , and let M' be the 
unique initial monomial ideal of the wall ideal which contains x v . If M' 
lies on HilbA then {M, M'} is an edge of the flip graph. We now illustrate 
the Macaulay 2 procedure for computing all flip neighbors of a monomial 
A-graded ideal. 

122 : f indPositiveVector = (m,s) -> ( 

expvector := first exponents s - first exponents m; 
n := #expvector; 

i := first positions(0. .n-1, j -> expvector_j > 0); 
splice {i:0, 1, (n-i-l):0> 

123 : flips = (M) -> ( 

R := ring M; 

— store generators of M in monoms 
monoms := first entries generators M; 
result := {}; 

— test each generator of M to see if it leads to a neighbo • • • 
scan(#monoms , i -> ( 

m : = monoms_i ; 

rest := drop (monoms , {i , i}) ; 

b := basis (degree m, R) ; 

s := (compress (b '/, M))_(0,0); 

J := ideal (m-s) + ideal rest; 

if poincare coker gens J == poincare coker gens M then ( 
w := f indPositiveVector(m,s) ; 

Rl := (coef f icientRing R) [generators R, Weights=>w] ; 

J = substituted, Rl) ; 

J = trim ideal leadTerm J; 

result = append (result, J) ; 

))); 

scan(#result , i->(print gens result_i))); 

The above code inputs a monomial ^4-graded ideal M whose minimal 
generators are stored in the list monoms. The flip neighbors of M will be 
stored in result. For each monomial x u in monoms we need to test whether 
it yields a flip neighbor of M or not. At the i-th step of this loop, we let 
m be the i-th monomial in monoms. The list rest contains all monomials in 
monoms except m. We compute the standard monomial s of M of the same 
degree as m. The wall ideal of m — s is the binomial ideal J generated by 
m — s and the monomials in rest. We then check whether J is ^4-graded by 
comparing its Hilbert series with that of M. If this is the case, we use the 
subroutine f indPositiveVector to find a unit vector w = (0, . . . , 1, . . . ,0) 
such that w ■ s > w ■ m. The flip neighbor is then the initial ideal of J with 
respect to w and it is stored in result. The program outputs the minimal 
generators of each flip neighbor. Here is an example. 
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124 : R = QQ [a. . e ,Degrees=>transpose A]; 

125 : M = ideal(a*e,c*d,a*c,a~2*d~2,a~2*b*d,a~3*d,c~2*e"3, 

c~3*e~2,c~4*e,c~5,c*e~5,a*d~5,b*e~6) ; 

o25 : Ideal of R 

126 : flips M 

{0} I ae cd ac a2d2 a3d c4 c2e3 c3e2 ad5 ce5 be6 I 

{0} I cd ae ac a2d2 a2bd a3d c3e2 c4e c5 ad4 ce5 c2e4 be6 I 

{0} I ae cd ac a2d2 a3d a2bd c2e3 c3e2 c4e c5 ce5 bce4 ad6 be6 I 

{0} I ae ac cd a2bd a3d a2d2 c2e3 c3e2 c4e c5 ce5 ad5 d7 I 

It is an open problem whether the toric Hilbert scheme HilbA is con- 
nected. Recent work in geometric combinatorics [|l7| suggests that this is 
probably false for some A. This result and its implications for HilbA will 
be discussed further in Section 2. However, we have the following result of 
Maclagan and Thomas |ll 

Theorem 1.2. The flip graph of A is connected if and only if the toric 
Hilbert scheme HilbA is connected. 

We now have two algorithms for listing monomial ideals on HilbA- First, 
there is the backtracking algorithm whose Macaulay 2 implementation was 
described above. Second, there is the flip search algorithm which starts with 
any coherent monomial ideal M and then constructs the connected compo- 
nent of M in the flip graph of A by carrying out local flips as above. This 



procedure is also implemented in TiGERS [10]. Clearly, the two algorithms 
will produce the same answer if and only if HilbA is connected. In other 
words, finding an example where HilbA is disconnected is equivalent to find- 
ing a matrix A for which the flip search algorithm produces fewer monomial 
ideals than the backtracking algorithm. 



2. Polyhedral Geometry 

Algorithms from polyhedral geometry are essential in the study of the 
toric Hilbert scheme. Consider the problem of deciding whether or not a 
given monomial ideal M in HilbA is coherent. This problem gives rise to a 
system of linear inequalities as follows. Let x Ul , . . . , x Ur be the minimal gen- 
erators of M, and let x Vi be the unique standard monomial with Aui = Av{. 
Then M is coherent if and only if there exists a vector w £ W 1 such that 
w ■ (ui — Vi) > for i = 1, . . . , r. Thus the test for coherence amounts 
to solving a feasibility problem of linear programming, and there are many 
highly efficient implementations (based on the simplex algorithms or inte- 
rior point methods) available for this task. For our experimental purposes, 
it is convenient to use the code polarCone .m2, given in Appendix 1, which 
is based on the (inefficient but easy-to-implement) Fourier- Motzkin elimi- 
nation method f23|| . This code converts the generator representation of a 
polyhedron to its inequality representation and vice versa. A simple exam- 
ple is given in Appendix 1. In particular, given a Grobner basis Q of I a, the 
function polarCone will compute all the extreme rays of the Grobner cone 
{w G R n : w ■ (ui -Vi)>0 for each x Ul - x Vi G Q}. 
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We now show how to use Macaulay 2 to decide whether a monomial A- 
graded ideal M is coherent. The first step in this calculation is to compute all 
the standard monomials of M of the same degree as the minimal generators 
of M. We do this using the procedure stdMonomials. 

127 : stdMonomials = (M) -> ( 

R := ring M; 
RM := R/M; 

apply (numgens M, i -> ( 

s := basis (degree (M_i) ,RM) ; lift(s_(0,0) , R))) 

) ! 

As an example, consider the following monomial A-graded ideal. 

128 : R = QQ [a. . e .Degrees => transpose A ]; 

129 : M = ideal(a"3*d, a~2*b*d, a~2*d~2, a*b~3*d, a*b"2*d"2, a*b*d~3, 

a*c, a*d~4, a*e, b~5*d, b"4*d"2, b~3*d~3, b~2*d~4, 
b*d"5, b*e, c*e"5) ; 

029 : Ideal of R 

130 : toString stdMonomials M 

030 = {b*c~3, c~4, c~3*e, c~5, c~4*e, c~3*e~2, b~2, c~2*e~3, b*d, c'6, • • • 
o30 : String 

From the pairs (x u , x v ) of minimal generators x u and corresponding stan- 
dard monomials x v , the function inequalities creates a matrix whose 
columns are the vectors u — v. 

131 : inequalities = (M) -> ( 

stds := stdMonomials (M) ; 

transpose matrix apply (numgens M, i -> ( 
flatten exponents (M_i) - 

flatten exponents (stds_i) ))) ; 

132 : inequalities M 

o32=|322 1 1 1 1 1 10000000 
-1 10321 -2 -1 5432110 
-3 -4 -3 -5 -4 -3 1 -2 -6 -5 -4 -3 -2 -1 1 
11212304 -1 12345 -1 -6 
-10 -1-2 -3 1 -1-2-3-4 1 5 

5 16 
o32 : Matrix ZZ < ZZ 

It is convenient to simplify the output of the next procedure using the 
following program to divide an integer vector by the g.c.d. of its components. 
We also load polarCone .m2 which is needed in decideCoherence below. 

133 : primitive := (L) -> ( 

n := #L-1; g := L#n; 

while n > do (n = n-1; g = gcd(g, L#n);); 
if g === 1 then L else apply(L, i -> i // g)); 

134 : load "polarCone .m2" 

135 : decideCoherence = (M) -> ( 

ineqs := inequalities M; 

c := first polarCone ineqs; 

m := - sum(numgens source c, i -> c_{i}) ; 

prods := (transpose m) * ineqs; 

if numgens source prods != numgens source compress prods 
then false else primitive (first entries transpose m)); 

Let K be the cone {x £ M. n : g ■ x < 0, for all columns g of ineqs }. The 
command polarCone ineqs computes a pair of matrices P and Q such that 
K is the sum of the cone generated by the columns of P and the subspace 
generated by the columns of Q. Let m be the negative of the sum of the 
columns of P. Then m lies in the cone —K. The entries in the matrix prods 
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are the dot products g-m for each column g of ineqs. Since M is a monomial 
A- graded ideal, it is coherent if and only if K is full dimensional which is 
the case if and only if no dot product g ■ m is zero. This is the conditional 
in the if . . then statement of decideCoherence. If M is coherent, the 
program outputs the primitive representative of m and otherwise returns the 
boolean false. Notice that if M is coherent, the cone —K is the Grobner 
cone corresponding to M and the vector m is a weight vector w such that 
in w (lA) = M. We now test whether the ideal M from line i29 is coherent. 

136 : decideCoherence (M) 
o36 = {0, 0, 1, 15, 18} 

036 : List 

Hence, M is coherent: it is the initial ideal with respect to the weight 
vector w = (0,0,1,15,18) of the toric ideal in our running example. The 
matrix in (jl|) has 55 noncoherent monomial ^4-graded ideals in total and 
here is one of them. 

137 : N = ideal(a*e,c*d,a*c,c~3*e,a~3*d,c~4,a*d~4,a~2*d~3,c*e~5, 

c~2*e-4,d-7) ; 

037 : Ideal of R 

138 : decideCoherence (N) 

038 = false 

In the rest of this section, we study the connection between A-graded 
ideals and polyhedral complexes defined on A which will relate the flip graph 
to the Baues graph of A. (See [jl6| for a survey of the Baues problem and 
its relatives). Let pos(A) := {Au : u £ M n ,u > 0} be the cone generated 
by the columns of A in M d . A polyhedral subdivision A of A is a collection 
of full dimensional subcones pos{A a ) of pos{A) such that the union of these 
subcones is pos(A) and the intersection of any two subcones is a face of 
each. Here A a := {a,- : j E a C {1, . . . ,n}}. It is customary to identify A 
with the set of sets {a : pos(A a ) £ A}. If every cone in the subdivision A 
is simplicial (the number of extreme rays of the cone equals the dimension 
of the cone), we say that A is a triangulation of A. The simplicial complex 
corresponding to a triangulation A is uniquely obtained by including in A 
all the subsets of every a £ A. We refer the reader to [20, §8] for more 
details. 

For each a E A, let I a be the prime ideal that is the sum of the toric ideal 
Ia„ and the monomial ideal (xj : j o~). Recall that (C*) n acts on R by 
scaling variables : A i— ► A • x := (Ai^i, . . . , X n x n ). Two ideals J and J' are 
said to be torus isomorphic if J = A • J' for some A € (C*) n . The following 
theorem shows that polyhedral subdivisions of A are related to A-graded 
ideals via their radicals. 



Theorem 2.1. (Theorem 10.10 |20], §10],) // / is any A-graded ideal, then 
there exists a polyhedral subdivision A(J) of A such that \fl = r\ ae /^mJ a 
where each component J a is a prime ideal that is torus isomorphic to I a . 
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We say that A(J) supports the A-graded ideal /. When M is a monomial 
^4-graded ideal, A(M) is a triangulation of A. In particular, if M is coherent 
(i.e, M = in w {lA) for some weight vector w), then A(M) is the regular or 
coherent triangulation of A induced by w [p0| , §8]. The coherent triangula- 
tions of A are in bijection with the vertices of the secondary polytope of A 

It is convenient to represent a triangulation A of A by its Stanley-Reisner 
ideal /a := {x^x^ ■ ■ ■ Xi k : {ii,«2j... , ifc} is a non-face of A). If M is a 



monomial ^4-graded ideal, Theorem 2.1 implies that I mm) ls the radical of 



M. Hence, we will represent triangulations of A by their Stanley-Reisner 
ideals. As seen below, the matrix in our running example has eight distinct 
triangulations corresponding to the eight distinct radicals of the 281 mono- 
mial A- graded ideals computed earlier. All eight are coherent. 

{{1,2}, {2, 3}, {3, 4}, {4, 5}} <-» (ac,ad,ae,bd,be,ce) 

{{1,3}, {3, 4}, {4, 5}} <-» (b,ad,ae,ce) 

{{1,2}, {2, 4}, {4, 5}} <- (c,ad,ae,be) 

{{1,2}, {2, 3}, {3, 5}} <-» (d,ac,ae,be) 

{{1,3}, {3, 5}} <-> (b,d,ae) 

{{1,4}, {4, 5}} <-> (b,c,ae) 

{{1,2}, {2, 5}} <- (c,d,ae) 

{{1,5}} <- (b,c,d) 

The Baues graph of A is a graph on all the triangulations of A in which 



two triangulations are adjacent if they differ by a single bistellar flip [16]. 
The Baues problem from discrete geometry asked whether the Baues graph 
of a point configuration can be disconnected for some A. Every edge of 
the secondary polytope of A corresponds to a bistellar flip and hence, the 
subgraph of the Baues graph that is induced by the coherent triangulations 
of A is indeed connected: it is precisely the edge graph of the secondary 
polytope of A. The Baues problem was recently settled by Santos |l7j who 
gave an example of a six dimensional point configuration with 324 points 
for which there is an isolated (necessarily non-regular) triangulation. 

Santos' configuration would also have a disconnected flip graph and hence 
a disconnected toric Hilbert scheme if it were true that every triangulation 
of A supports a monomial A-graded ideal. However, Peeva has shown that 
this need not be the case (Theorem 10.13 in [^, §10]). Hence, the map from 
the set of all monomial A-graded ideals to the set of all triangulations of A 
that sends M i— > A(M) is not always surjective, and it is unknown whether 
Santos' 6 x 324 configuration has a disconnected toric Hilbert scheme. 

Thus, even though one cannot in general conclude that the existence of a 
disconnected Baues graph implies the existence of a disconnected flip graph, 
there is an important special situation in which such a conclusion is possible. 
We call an integer matrix A of full row rank unimodular if the absolute value 
of each of its non-zero maximal minors is the same constant. A matrix A is 
unimodular if and only if every monomial j4-graded ideal is square-free. For 
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a unimodular matrix A, the Baues graph of A coincides with the flip graph 
of A, and as you might expect, Santos' configuration is not unimodular. 

Theorem 2.2. (Lemma 10. 14 |2(], §10] j If A is unimodular, then each tri- 
angulation of A supports a unique (square-free) monomial A-graded ideal. 
Moreover, a monomial A-graded ideal is coherent if and only if the triangu- 
lation supporting it is coherent. 
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Theorem 2.2 provides a competitive algebraic algorithm for computing 
all the triangulations of a unimodular matrix since they are precisely the 
polyhedral complexes supporting monomial A-graded ideals. Then we could 
enumerate the connected component of a coherent monomial A-graded ideal 
in the flip graph of A to decide whether the Baues/flip graph is disconnected. 

Let A r be the standard r-simplex which is the convex hull of the r + 1 
unit vectors in W +1 and let A(r, s) be the (r + s + 2) x (r + 1) (s + 1) matrix 
whose columns are the products of the vertices of A r and A s . All matrices 
of type A(r, s) are unimodular. From the product of two triangles we get 



4(2,2) 



We can now use our algebraic algorithms to compute all the triangulations 
of A(2, 2). Since Macau/ay 2 requires the first entry of the degree of every 
variable in a ring to be positive, we use the following matrix with the same 
row space as A(2, 2) for our computation: 

139 : A22 = 

{{1,1, 1,1, 1,1, 1,1,1} ,{0,0,0, 1,1, 1,0, 0,0}, {0,0, 0,0, 0,0, 1,1,1}, 
{1,0, 0,1, 0,0, 1,0,0}, {0,1, 0,0, 1,0, 0,1,0}, {0,0, 1,0, 0,1, 0,0,1}}; 

140 : 122 = toricldeal A22 

o40 = ideal (f*h - e*i, c*h - b*i, f*g - d*i, e*g - d*h, c*g - a*i, b* • • • 

040 : Ideal of R 

The ideal 122 is generated by the 2 by 2 minors of a 3 by 3 matrix of 
indeterminates. This is the ideal of P 2 x P 2 embedded in P 8 via the Segre 
embedding. 

141 : Graver22 = graver 122; 

1 15 

041 : Matrix R < R 

142 : generateAmonos (Graver22) ; 

143 : #monos 
o43 = 108 

144 : scan(0 . . 9, i->print toString monos#i) 
{f*h, c*h, f*g, e*g, c*g, b*g, c*e, c*d, b*d} 
{f*h, d*h, c*h, f*g, c*g, b*g, c*e, c*d, b*d} 
{d*i, f*h, d*h, c*h, c*g, b*g, c*e, c*d, b*d} 
{e*i, c*h, f*g, e*g, c*g, b*g, c*e, c*d, b*d} 
{e*i, d*i, c*h, e*g, c*g, b*g, c*e, c*d, b*d} 
{e*i, d*i, d*h, c*h, c*g, b*g, c*e, c*d, b*d} 
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{f*h, c*h, f*g, e*g, c*g, b*g, c*e, a*e, c*d} 
{e*i, c*h, f*g, e*g, c*g, b*g, c*e, a*e, c*d, b*d*i} 
{e*i, c*h, f*g, e*g, c*g, b*g, c*e, a*e, c*d, a*f*h} 
{e*i, d*i, c*h, e*g, c*g, b*g, c*e, a*e, c*d} 

Thus there are 108 monomial A{2, 2)-graded ideals and decideCoherence 
will check that all of them are coherent. Since A(2, 2) is unimodular, each 
monomial A{2, 2)-graded ideal is square-free and is hence radical. These 
108 ideals represent the 108 triangulations of A(2, 2) and we have listed 
ten of them above. The flip graph (equivalently, Baues graph) of A{2, 2) is 
connected. However, it is unknown whether the Baues graph of A(r, s) is 
connected for all values of (r, s). 

3. Local Equations 
Consider the reduced Grobner basis of a toric ideal I a for a term order 

w: 

(2) { x Ul - x Vl , x U2 -x v \ ... , x Ur -x Vr }. 

The initial ideal M = in w {lA) = (x ui , x U2 , . . . , x Ur ) is a coherent monomial 
^4-graded ideal. In particular, it is a (C*) n -fixed point on the toric Hilbert 



scheme HilbA- We shall explain a method, due to Peeva and Stillman [14], 
for computing local equations of HilbA around such a fixed point. A variant 
of this method also works for computing the local equations around a non- 
coherent monomial ideal M, but that variant involves local algebra, specifi- 
cally Mora's tangent cone algorithm, which is not yet fully implemented in 



Macaulay 2 . See [14] for details. 

We saw how to compute the flip graph of A in Section 1. The vertices 
of this graph are the (C*) Tl -fixed points M and its edges correspond to the 
(C*) n -fixed curves. By computing and decomposing the local equations 
around each M, we get a complete description of the scheme HilbA- 

The first step is to introduce a new variable Zi for each binomial in our 
Grobner basis (g) and to consider the following r binomials: 

(3) X U1 - Z X ■ X V1 , X U2 - Z 2 ■ X V2 , . . . , X Ur - Z r ■ X Vr 

in the polynomial ring C[x,z] in n + r indeterminates. The term order w 
can be extended to an elimination term order in C[x,z] so that x Ui is the 
leading term of x Ui — Z{ ■ x Vi for all i. We compute the minimal first syzygies 
of the monomial ideal M, and form the corresponding S'-pairs of binomials 
in (||). For each S'-pair 

lcm(x u %x u A . ... ... lcm(x u \x u i) , ... .... 

■ • (x Ul - Zi-X" 1 ) - ' • (X^ - Zj ■ x V] ) 

we compute a normal form with respect to (Q) using the extended term order 
w. The result is a binomial in C[x,z] which factors as 

af*.zP. {zi - z s ), 

where a £ N n and /3, 7,5 G N r . Note that this normal form is not unique 
but depends on our choice of a reduction path. Let Jm denote the ideal in 
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C[zi, ... ,z r ] generated by all binomials z 13 ■ (z 7 — z s ) gotten from normal 
forms of all the S'-pairs considered above. 

Proposition 3.1. (Peeva-Stillman [|14|]) The ideal Jm is independent of 
the reduction paths chosen. It defines a subscheme of C r isomorphic to an 
affine open neighborhood of the point M on the toric Hilbert scheme HilbA- 

We apply this technique to compute a particularly interesting affine chart 
of HilbA for our running example. Consider the following set of 13 binomials: 

{ ae — z\bd, cd — Z2be, ac — z^b 2 , a 2 d 2 — z^c 3 e, a 2 bd — z§c 4 , 

a 3 d — z 6 bc 3 , c 2 e 3 — z 7 ad 4 , c 3 e 2 — z 8 abd 3 , c 4 e — zgab 2 d 2 , 

c 5 — zioab 3 d, ce 5 — znd 6 , ad 5 — z^bce 1 , be G — z^d 7 }. 

If we set z\ = Z2 = ■ ■ ■ = z\z = 1 then we get a generating set for the toric 
ideal I a- The 13 monomials obtained by setting z\ = z-i = ■ ■ ■ = z±s = 
generate the initial monomial ideal M = in w {lA) with respect to the weight 
vector w = (9,3,5,0,0). Thus M is one of the 226 coherent monomial A- 
graded ideals of our running example. The above set of 13 binomials in 
C[x, z] give the universal family for HilbA around this M. 

The local chart of HilbA around the point M is a subscheme of affine space 
C 13 with coordinates Z\, . . . , Z13, whose defining equations are obtained as 
follows. Extend the weight vector w by assigning weight zero to all variables 
Zi, so that the first term in each of the above 13 binomials is the leading 
term. For each pair of binomials corresponding to a minimal syzygy of M, 
form their S-pair and then reduce it to normal form with respect to the 13 
binomials above. For instance, 

£(c 5 — z±oab 3 d, ce 5 — znd 6 ) = znc 4 d 6 —zioab s de 5 — > b 4 d 2 e 4 -(z2Zu—ziZio). 

Each such normal form is a monomial in a, b, c, d, e times a binomial in 
Z\ , . . . , Z13. The set of all these binomials generates the ideal Jm of local 
equations of HilbA around M. In our example, Jm is generated by 27 
nonzero binomials. This computation can be done in Macaulay 2 using the 
procedure localCoherentEquations. 

i45 : localCoherentEquations = (IA) -> ( 

— IA is the toric ideal of A living in a ring equipped 

— with weight order w, if we are computing the local 

— equations about the initial ideal of IA w.r.t. w. 
R := ring IA; 

w := (monoid R) .Options. Weights; 

M := ideal leadTerm IA; 

S := first entries ((gens M) '/, IA) ; 

— Make the universal family J in a new ring, 
nv : = numgens R ; n : = numgens M ; 

T = (coef f icientRing R) [generators R, z_l . . z_n, 

Weights => flatten splice{w, n:0}, 
MonomialSize => 16] ; 

M = substitute (generators M,T); 

S = apply(S, s -> substituted, T)) ; 

J = ideal apply (n, i -> 

M_(0,i) - T_(nv + i) * S_i) ; 

— Find the ideal Ihilb of local equations about M: 
spairs := (gens J) * (syz M) ; 

g := forceGB gens J; 

B = (coef f icientRing R) [z_l . . z_n,MonomialSize=>16] ; 
Fones := map(B,T, matrix(B , {splice {nv:l}}) I vars B) ; 
Ihilb := ideal Fones (spairs g) ; 
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Ihilb 

); 

Suppose we wish to calculate the local equations about M = in w {lA)- 
The input to localCoherentEquations is the toric ideal Ia living in a 
polynomial ring equipped with the weight order specified by w. This is 
done as follows. 

146 : IA = toricldeal A; 
o46 : Ideal of R 

147 : Y = QQ[a..e, MonomialSize => 16, 

Degrees => transpose A, Weights => {9,3,5,0,0}]; 

148 : IA = substituted^, Y); 

048 : Ideal of Y 

The initial ideal M is calculated in the third line of the algorithm, and 
S stores the standard monomials of M of the same degrees as the mini- 
mal generators of M. We could have calculated S using our old procedure 
stdMonomials but this involves computing the monomials in i?^ for var- 
ious values of b which can be slow on large examples. As by-products, 
localCoherentEquations also gets J, the ideal of the universal family for 
HilbA about M, the ring T of this ideal, and the ring B of Ihilb which is 
the ideal of the affine patch of HilbA about M. The matrix spairs contains 
all the S'-pairs between generators of J corresponding to the minimal first 
syzygies of M. The command f orceGB is used to declare the generators of 
J to be a Grobner basis, and Fones is the ring map from T to B that sends 
each of a, b, c, d, e to one and the z variables to themselves. The columns 
of the matrix (spairs % g) are the normal forms of the polynomials in 
spairs with respect to the forced Grobner basis g and the ideal Ihilb of 
local equations is generated by the image of these normal forms in the ring 
B under the map Fones. 

149 : JM = localCoherentEquations (IA) 

049 = ideal (zz -z,zz -z,-zz +z,-zz +z,-zz + ••• 

12 312 3 47 2 58 2 15 

o49 : Ideal of B 

Removing duplications among the generators, 

Jm = (zi - z w zu,z 2 - z A z 7 , z 2 - z 5 z 8 , z 2 - ZuZ 12 , z 2 - Z 1 Z 11 Z 13 , 

Z3 - Z±Z 2 , Z 3 - Z 5 Z 9 , Z 4 - ZlZ 5 , Zq - Z 3 Z 5 , Zq - ZlZ 2 Z 5 , z 7 - Z 1 Z W , z 8 - z 1 z 7 , 
Zc, - ZiZs, Z12 - ZiZis, Z 1 Z 2 - Z 5 Zg, Z\Z 2 - ZiZ 5 Z$, Z\Z 2 - zjz4Zi , Z\Z 2 - zfz 5 Z 7 , 

z\z 2 - zizuz 12 , z 1 z 2 - z 2 z w zn,zfz 4: - z 3 z 11 ,z 1 z 5 z 8 - z 4 z 8 , z 2 z w - ziz i2 , 

Z3Z4 - ZlZQ, Z 3 Z 7 - Z 2 Z 8 , Z 3 Z$ - Z 2 Z 9 , Z3Z10 - z 2 z 7 ). 

Notice that there are many generators of Jm which have a single vari- 
able as one of its terms. Using these generators we can remove variables 
from other binomials. This is done in Macaulay 2 using the subroutine 
removeRedundantVariables which is the main ingredient of the package 
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minPres . m2 for computing the minimal presentations of polynomial quotient 
rings. Both removeRedundant Variables and minPres.m2 are explained 
in Appendix 2. The command removeRedundantVariables applied to an 
ideal in a polynomial ring (not quotient ring) creates a ring map from the 
ring to itself that sends the redundant variables to polynomials in the non- 
redundant variables and the non-redundant variables to themselves. Apply- 
ing this to our ideal Jm we obtain the following simplifications. 

150 : load "minPres .m2" ; 

151 : G = removeRedundantVariables JM 

32 43 2 4 3 2 • • • 

o51 = map(B,B,-[z z , z z z , z z z , z z z , z , z z z , z ••• 
10 11 5 10 11 5 10 11 5 10 11 5 5 10 11 10 • • • 

051 : RingMap B < B 

152 : ideal gens gb(G JM) 

3 2 2 

052 = ideal (z z z - z z z ) 

5 10 11 10 11 13 

052 : Ideal of B 

Thus our affine patch of HilbA has the coordinate ring 

^ r 1/T C[z 5 ,z w ,z u ,z 13 ] C[z 5 ,z w ,z n ,z 13 ] 

<L [z\,Z2, ■ ■ ■ ,213 \/Jm — -, 5 — o 2 r = 77 o s 2~T- 

\ z $ z io z ii - zioznzn) {{z 5 zf - 213)^10%) 

Hence, we see immediately that there are three components through the 
point M on HilbA- The restriction of the coherent component to the 
affine neighborhood of M on HilbA is defined by the ideal quotient (Jm ■ 
(z\Z2 ■ ■ ■ Z13) 00 ) and hence the first of the above components is an affine 
patch of the coherent component. Locally near M it is given by the single 
equation z^zf — z± 3 = in A 4 . It is smooth and as expected, has dimension 
three. The second component, z±q = is also of dimension three and is 
smooth at M. The third component, given by z\^ = is more interesting. 
It has dimension three as well, but is not reduced. Thus we have proved the 
following result. 

Proposition 3.2. The toric Hilbert scheme HilbA of the matrix A in ([Jj 
is not reduced. 

We can use the ring map G from above to simplify J so as to involve only 
the four variables 25, z\q, Z\\ and Z\ 3 . 

153 : CX = G)Q[a. . e, z_5,z_10,z ll,z_13, Weights => splice 

{9,3,5,0,0,0,0,0,0}]; 

154 : F = map(CX, ring J, matrix{{a,b,c,d,e}} I 

substitute (G .matrix, CX) ) 

3 2 4 3 

o54 = map(CX,T,{a, b, c, d, e, z z ,zz z ,zz z ,zz z , z • • • 

10 11 5 10 11 5 10 11 5 10 11 

054 : RingMap CX < T 

Applying this map to J we get the ideal Jl, 

155 : Jl = F J 
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3 2 2 4 3 

o55 = ideal (c*d - b*e*z z , a*e - b*d*z z z , a*c - b z z z , a • • • 
10 11 5 10 11 5 10 11 

055 : Ideal of CX 

and adding the ideal (zf-^) to Jl we obtain the universal family for the 
non-reduced component of HilbA about M. 

i56 : substitute(ideal(z_ll~2) ,CX) + Jl 

2 3 2 2 4 

056 = ideal (z , c*d - b*e*z z , a*e - b*d*z z z , a*c - b z z z • • • 

11 10 11 5 10 11 5 10 • • • 

056 : Ideal of CX 

In the rest of this section, we present an interpretation of the ideal Jm in 
terms of the combinatorial theory of integer programming. See, for instance, 



0, §4] or [22 1 for the relevant background. Our reduced Grobner basis 



is the minimal test set for the family of integer programs 

(4) Minimize w ■ u subject to A ■ u = b and u G N n , 

where A E N dxn and w G Z n are fixed and b ranges over N d . If u' G N ra is 
any feasible solution to , then the corresponding optimal solution u G N n 
is computed as follows: the monomial x u is the unique normal form of x u 
modulo the Grobner basis (g). 

Suppose we had reduced x u modulo the binomials (H|) instead of (||). 
Then the output has a z-factor which depends on our choice of reduction 
path. To be precise, suppose the reduction path has length m and at the 
j-th step we had used the reduction x ^ — ► z^. ■ x" 11 ^ . Then we would 
obtain the normal form 



Reduction paths can have different lengths. If we take another path which 
has length m! and uses x Uv i — > z v . ■ x Vv i at the j'-th step, then the output 
would be 

Z V\ Z V2 Z V3 ' ' ' Z V m l ' X ■ 

Theorem 3.3. The ideal Jm of local equations on HilbA is generated by 
the binomials 

Z H1 Z ^2 Z ^3 ' ' ' Z ^m ~ Z V\ Z V2 Z V3 ' ' ' Z V m l 

each encoding a pair of distinct reduction sequences from a feasible solution 
of (j^) to the corresponding optimal solution using the minimal test set in 

Proof: The given ideal is contained in Jm because its generators are differ- 
ences of monomials arising from the possible reduction paths of lcm(x Ui , x U: > '), 
for 1 < i, j < r. Conversely, any reduction sequence can be transformed into 
an equivalent reduction sequence using S-pair reductions. This follows from 
standard arguments in the proof of Buchberger's criterion 0, §2.6, Theorem 
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6], and it implies that the binomials ■ ■ ■ z^ m — z Ul - - - z v , are C[z] -linear 
combinations of the generators of Jm ■ D 

A given feasible solution of the integer program (||) usually has many 
different reduction paths to the optimal solution using (0). For our matrix 
(ID) and cost vector w = (9,3,5,0,0), the monomial a bde 6 encodes the 
feasible solution (2,1,0,1,0,6) of the integer program (Q) with right hand 
side vector ( 5g ). There are 19 different paths from this feasible solution 
to the optimal solution (0,3,0,3,4) encoded by the monomial 6 3 d 3 e 4 . The 
generating function for these paths is: 

z\ + 3z 1 z 2 z 5 z 7 + 2z 1 z 2 z 5 z 2 ,z 12 + 2z 1 z 2 z 5 z 8 + 2z 1 z 2 z 12 z 13 + z 1 z 5 z 9 

+ Z%Z4,Z 5 Z% + z\z±Z 13 + z%z 5 Zu + 2z 2 z 3 z 5 z 7 + Z 3 Z 5 Z^Z 12 + Z 3 Z 5 Z 8 + Z 3 Z 12 Z 13 . 

The difference of any two monomials in this generating function is a valid 
local equation for the toric Hilbert scheme of ([j]). For instance, the binomial 
z 3 Z5Zj z\ 2 — z 3 zi 2 zi 3 lies in Jm, and, conversely, Jm is generated by binomials 
obtained in this manner. 

The scheme structure of Jm encodes obstructions to making certain re- 
ductions when solving our family of integer programs. For instance, the 
variable z 3 is a zero-divisor modulo Jm- If we factor it out from the bino- 
mial z 3 zz,ZjZi 2 — z 3 z\ 2 z\ 3 6 Jm , we get zz,z 7 zi 2 — z\ 2 z\ 3 which does not lie 
in Jm- Thus there is no monomial a ll b l2 c l3 d u e l5 for which both the paths 
Z5Z 7 z\ 2 and z\ 2 z\ 3 are used to reach the optimum. It would be a worthwhile 
combinatorial project to study the path generating functions and their re- 
lation to the ideal Jm in more detail. 

It is instructive to note that the binomials z^ z^ 2 - - - z^ m — z Ul z V2 ■ ■ ■ z v , 
in Theorem [Di] do not form a vector space basis for the ideal Jm- We 
demonstrate this for the lexicographic Grobner basis (with a y- b >- c >- 
d y- e) of the toric ideal defining the rational normal curve of degree 4. In 

this case, we can take A = ( J \ \ J \ ) and the universal family in 
question is : 

C O OO O 1 

\ac — z\b , ad — z 2 bc, ae — z 3 c , bd — z^c , be — z$cd, ce — z%d ). 

The corresponding ideal of local equations is Jm = (-23 — ^2-25, z 2 — Z\Z±, z$ — 
Z4ZQ), from which we see that M is a smooth point of HiWa- The binomial 
z\z<z — Z1Z4ZQ lies in Jm but there is no monomial which has the reduction 
path z\z§ or z^z\ to optimality. Indeed, any monomial which admits the 
reductions z\z^ or z^z\ must be divisible by either ace or abe. The path 
generating functions for these two monomials are 

abe — > (z 3 + ziz 4 z 5 + z 2 z 5 ) ■ be 2 

ace — > (z 3 + z\ Z4Z5 + z 2 Z4Zq) • c 3 . 

Thus every reduction to optimality using z\ and 25 must also use Z4, and 
we conclude that z\z§ — z\z^z% is not in the C-span of the binomials listed 



in Theorem 3.2 
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4. The coherent component of the toric Hilbert scheme 

In this section we study the distinguished component of the toric Hilbert 
scheme HilbA- The toric ideal Ia is a point on this component. We will 
show that this component need not be normal, and will describe how its 
local and global equations can be computed using Macaulay 2 . Every term 
order for the toric ideal Ia can be realized by a weight vector which is an 
element in the lattice N = Homz(kerz(A),lj) ~ 7L n ~ d . Two weight vectors 
w and w' in N are considered equivalent if they define the same initial ideal 
in w (lA) = in w i(lA)- These equivalence classes are the relatively open cones 



of a projective fan in the lattice N called the Grobner fan of Ia [13], 



[211- 



Theorem 4.1. The toric ideal I a lies on a unique irreducible component 
of the toric Hilbert scheme HilbA, called the coherent component. The nor- 
malization of the coherent component is the projective toric variety defined 
by the Grobner fan of I a ■ 

Proof: The divisor at infinity on the toric Hilbert scheme HilbA consists 
of all points at which at least one of the local coordinates (around some 
monomial A-graded ideal) is zero. This is a proper Zariski closed subset of 
HilbA, consisting of all those A-graded ideals which contain at least one 
monomial. By Lemma 10.12], its complement in HilbA consists of 
precisely the orbit of I a under the action of the torus (C*) n . This shows that 
the closure of the (C*) n -orbit of I a is an irreducible component of HilbA- 
We call it the coherent component. 

Identifying (C*) n with ffomz(Z n , C*), we note that the stabilizer of Ia 
consists of those linear forms w which restrict to zero on the kernel of A. 
Therefore the coherent component is the closure in HilbA of the orbit of 
the point Ia under the action of the torus N (g> C* = Homz(kerz(A), C*). 
The (iV (g) C*)-fixed points on this component are precisely the coherent 
monomial ^4-graded ideals, and the same holds for the toric variety of the 
Grobner fan. 

Fix a maximal cone a in the Grobner fan S^, and let M = (x Ul , . . . , x Ur ) 
be the corresponding (monomial) initial ideal of I a- As before we write 

{X U1 - Zl ■ X Vl , X U * -Z 2 -X V \ ... , X Ur - Z r ■ X Vr ) 

for the universal family arising from the corresponding reduced Grobner 
basis of I a ■ Let Jm be the ideal in C [z\ , z 2 , ■ ■ • , z r ] defining this family. 

The restriction of the coherent component to the affine neighborhood 
of M on HilbA is defined by Jm '■ {z\z 2 ■ ■ ■ z r )°°. It then follows from 
our combinatorial description of the ideal Jm that this ideal quotient is a 
binomial prime ideal. In fact, it is the ideal of algebraic relations among the 
Laurent monomials x Ul ~ Vl , . . . ,x Ur ~ Vr . We conclude that the restriction of 
the coherent component to the affine neighborhood of M on HilbA equals 

(5) SpecC[x Ul - 
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The abelian group generated by the vectors u% — v\, . . . ,u r — v r equals 
kerz(A) = Homz(N, Z). The cone generated by the vectors Ui—vi,... ,u r — 
v r is precisely the polar dual c v to the Grobner cone a. Therefore the 
normalization of the affine variety (|5|) is the normal afflne toric variety 

(6) Spec C[ker z (A) n a v ] . 

The various normalization maps from @ to (|5|), obtained as a varies 
over the cones of £a are compatible with gluing along S^. The result is 
the desired normalization map from the projective toric variety associated 
with the Grobner fan of I a onto the coherent component of the toric Hilbert 
scheme HilbA- D 

We now present an example which shows that the coherent component 
of HilbA need not be normal. This example is derived from the matrix 
that appears in Example 3.15 of §. This example is also mentioned in 15] 
without details. Let d = 4 and n = 7 and fix the matrix 



/ 1 


1 


1 


1 


1 


1 


1 \ 





6 


7 


5 


8 


4 


3 


3 


7 


2 





7 


6 


1 


V 6 


5 


2 


6 


5 








The lattice N = Homz(kerz(A),Z) is three-dimensional. The toric ideal 
I A is minimally generated by 30 binomials of total degree between 6 and 93. 

157 : A = {{1,1, 1,1, 1,1,1}, {0,6, 7, 5, 8, 4, 3}, {3, 7, 2, 0,7, 6,1}, 

{6,5,2,6,5,0,0}}; 

158 : IA = toricldeal A 

23 32 2 44 84 436 724 4... 

o58 = ideal (ace- b*d f , a c*d*e f -bg,def -beg, a*b c • • • 

058 : Ideal of R 

We fix the weight vector w = (0, 0, 276, 220, 0, 0, 215) in N and compute 
the initial ideal M = in w {lA)- This initial ideal has 44 minimal generators. 

159 : Y = QQ[a..g, MonomlalSlze => 16, 

Weights => {0,0,276,220,0,0,215}, 
Degrees =>transpose A] ; 

160 : IA = substitute(IA,Y) ; 

060 : Ideal of Y 

161 : M = ideal leadTerm IA 

23 84 724 473 5435 2654 3 3 1 • • • 

061 = ideal (ace,bg,bcg, a*b cf,bcdf,abcg,abc ••• 

o61 : Ideal of Y 

Proposition 4.2. The three dimensional affine variety for the initial 
ideal M with respect to w = (0, 0, 276, 220, 0, 0, 215) of the toric ideal of A 
in Q), is not normal. 
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Proof. The universal family for the toric Hilbert scheme HiWa at M is: 

{ a 2 e 15 5 18 _ Zlb 3 c 6 d 10 / 16 ) 6 13 d l 5/ 16 _ ^S^U ^STjllO _ ^92^ 

ac lA d ll f 23 -z^g 29 , b 7 c 2 g A - z b d A e 3 f\ bc 3A d 32 f 2 - z Ai e> 3 g™}. 

The semigroup algebra in (||) is generated by 44 Laurent monomials gotten 
from this family. It turns out that the first four monomials suffice to gen- 
erate the semigroup. In other words, for all j £ {5, 6, . . . , 44} there exist 
*1j*2j*3)*4 £ N such that Zj — z\ x z l 2 2 z l 3 3 z % £ £ Jm ■ (z± • • -244)°°. Hence the 
semigroup algebra in (||) is: 

r a 2 e 15 g 18 b 13 d 15 f 16 c 59 d 57 f no ac 1A d ll f 23 , C[ Zl , z 2 , z 3 , z 4 ] 

L ^3 c 6 d l 0/ 16' a 8 ce 21 5 14' e 92 5 134 ' h( AV g 29 J ~ ( Z \z 2 Z 3 - zf) ' 

This algebra is not integrally closed, since a toric hypersurface is normal if 
and only if at least one of the two monomials in the defining equation is 
square-free. Its integral closure in C[ker%(A)] is generated by the Laurent 
monomial 



z A . .1 b 5 c 26 d 3l f 



55 



2 



(8) 4 = {Z1Z2Z3 >~ ~ a^g<* 

Hence the affine chart (^) of the toric variety of the Grobner fan of I a is the 
spectrum of the normal domain C[zi, z 2 , z 3 ,y]/ (ziz 2 z 3 — y 2 ), where y maps 
to®. ^ 

We now examine the local equations of HilbA about M for this example. 

162 : JM = localCoherentEquations(IA) 

o62 = ideal (zz - z , z z -z,zz -z,zz - z , z z - z , z ••• 
12 3 12 3 15 4 13 6 13 6 1 • • • 

062 : Ideal of B 

163 : G = removeRedundantVariables JM; 

063 : RingMap B < B 

164 : toString ideal gens gb(G JM) 

064 = ideal(z_32*z_42~2*z_44-z_37~2*z_42,z_32~3*z_35*z_37~2-z_42~2*z_4 • • • 
o64 : String 

This ideal has six generators and decomposing it we see that there are 
five components through the monomial ideal M on this toric Hilbert scheme. 
They are defined by the ideals: 



(^32^42^44 — £37, ^32^35 ~ ^42, ^32^35^37 ~ 2 42 ^44, ^32 z 35^37 



^32^35^37 — zf 2 zl 4 , Z3^Z 37 — z\ 2 z\^) 



z 42 z 44' 



• (244, ^37 > 

• (^37,^42) 

• (^42,235) 

• (^42, ^| 2 )- 

All five components are three dimensional. The first component is an affine 
patch of the coherent component and two of the components are not reduced. 
Let K be the first of these ideals. 
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165 : K = ideal(z_32*z_42*z_44-z_37~2,z_32~4*z_35-z_42, 

z_32-3*z_35*z_37-2-z_42-2*z_44,z_32~2*z_35*z_37-4-z_42~3*z_4 • • • 
z_32*z_35*z_37~6-z_42~4*z_44-3,z_35*z_37-8-z_42-5*z_44~4) ; 

065 : Ideal of B 

Applying removeRedundantVariables to K we see that the affine patch 
of the coherent component is, locally at M, a non-normal hypersurface sin- 
gularity (agreeing with (|8|)). The labels on the variables depend on the order 
of elements in the initial ideal M computed by Macaulay 2 in line i61. 

166 : GG = removeRedundantVariables K; 

066 : RingMap B < B 

167 : ideal gens gb (GG K) 

5 2 

067 = ideal (z z z - z ) 

32 35 44 37 

o67 : Ideal of B 

There is a general algorithm due to de Jong || for computing the normal- 
ization of any affine variety. In the toric case, the problem of normalization 
amounts to computing the minimal Hilbert basis of a given convex ratio- 
nal polyhedral cone |18[| . An efficient implementation can be found in the 
software package Normaliz by Bruns and Koch Q. 

Our computational study of the toric Hilbert scheme in this chapter was 
based on local equations rather than global equations (arising from a pro- 
jective embedding of HilbA), because the latter system of equations tends 
to be too large for most purposes. Nonetheless, they are interesting. In the 
remainder of this section, we present a canonical projective embedding of 
the coherent component of HilbA- 

Let G\, G2, G3, • • • , G s denote all the Graver fibers of the matrix A. In 
Section 1 we saw how to compute them in Macaulay 2 . Each set Gi con- 
sists of the monomials in C[x] which have a fixed Graver degree. Consider 
the set G := G1G2G3 ■ ■ ■ G s which consists of all monomials which are 
products of monomials, one from each of the distinct Graver fibers. Let t 
denote the cardinality of G. We introduce an extra indeterminate z, and 
we consider the N-graded semigroup algebra C[zG], which is a subalgebra 
of C[x±, . . . ,x n ,z]. Labeling the elements of G with indeterminates yi, We 
can write 

C[zG] = C[ yi ,y 2 ,... ,y t ]/P A , 

where Pa is a homogeneous toric ideal associated with a configuration of t 
vectors in Z n+1 . We note that the torus (C*) n acts naturally on C[zG]. 

/ 3 2 1 \ 

Example 4.3. Let n = 4, d = 2 and A = I g)> so that I a is 

the ideal of the twisted cubic curve. There are five Graver fibers: 

168 : A = {{1,1, 1,1}, {0,1, 2, 3}}; 

169 : I = toricldeal A; 
o69 : Ideal of R 
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i70 



Graver = graver I ; 



o70 
i71 
i72 



1 5 
Matrix R < R 

fibers = graverFibers Graver; 

peek fibers 



o72 = MutableHashTable{{2, 2} => {0, 0} 

{2, 3} => {0, 0} 

{2, 4} => {0, 0} 

{3, 3} => {0, 0} 

{3, 6} => {0, 0} 



ac b2 I 
ad be I 
bd c2 I 
a2d abc b3 I 
ad2 bed c3 I 



072 : Net 

The set G = G1G2G3G4G5 consists of 22 monomials of degree 14. 

173 : G = trim product (values fibers, ideal) 

5 5 435 534 4224 34 4 264 4 • • ■ 

073 = ideal (a b*c*d ,abd,acd,abcd,ab c*d , a b d , a b • • • 

073 : Ideal of R 

174 : numgens G 

074 = 22 

We introduce a polynomial ring in 22 variables yi,V2,--- ,2/22) and we 
compute the ideal P4. It is generated by 180 binomial quadrics. 



i75 
i76 
i77 

o77 
i78 
i79 



z = symbol z ; 

S = QQ[a,b,c,d,z] ; 

zG = z ** substitute (gens G, S) ; 

1 22 
Matrix S < S 

R = QQ[y_l . . y_22] ; 

F = map(S,R,zG) 



5 5 435 534 4224 34 4 26 

o79 = map(S,R,{a b*c*d z, abdz, acdz, abedz, ab c*d z, a b 

o79 : RingMap S < R 

i80 : PA = trim ker F 



o80 = ideal (y - y y , y y -yy.yy -yy.yy 

21 20 22 19 21 18 22 18 21 17 22 17 21 

080 : Ideal of R 

These equations define a toric surface of degree 30 in projective 21-space. 

181 : codim PA 

081 = 19 

182 : degree PA 

082 = 30 

The surface is smooth, but there are too many equations and the codi- 
mension is too large to use the Jacobian criterion for smoothness ||, §16.6] 
directly. Instead we check smoothness for each open set yi ^ 0. 

183 : Aff = apply(1..22, v -> ( 

K = substitute (PA, y_v => 1); 

FF = removeRedundant Variables K; 

ideal gens gb (FF K))) ; 
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i84 : scan(Aff, i -> print toString i) ; 
ideal () 
ideal () 
ideal () 

ideal (y_l"4*y_5*y_21-l) 
ideal (y_l~4*y_6~6*y_21-l) 
ideal O 

ideal (y_l _ 2*y_ll-2*y_17-l) 
ideal (y_l"3*y_9-2*y_21-2-l) 

ideal (y_6"3*y_21-y_10,y_l*y_10"3-y_6"2,y_l*y_6*y_10"2*y_21-l) 
ideal (y_6*y_15-l,y_2*y_15~2-y_6*y_14,y_6"2*y_14-y_2*y_15) 
idealO 

ideal (y_ll*y_13-l,y_l _ 2*y_21-3-y_13-2) 
ideal (y_l"2*y_14~3*y_21 "3-1) 
ideal (y_10~ 2*y_21-l ,y_l*y_15~4-y_lCT 3) 
ideal 

ideal (y_ll*y_20-l,y_3*y_20~2-y_ll*y_17,y_ll~2*y_17-y_3*y_20) 
ideal(y_ll*y_18*y_21-l,y_l*y_21~3-y_ll*y_18~2,y_ll~2*y_18~3-y_l*y_21~2 • • • 
ideal (y_l*y_19-4*y_21-4-l) 
ideal(y_15*y_22-l) 
ideal 

ideal(y_20*y_22-l) 
ideal 

By examining these local equations, we see that HilbA is smooth, and 
also that there are eight fixed points under the action of the 2-dimensional 
torus. They correspond to the variables y%, y 2 , 2/3, 2/6, 2/11, 2/15, 2/20 and 2/22- 
By setting any of these eight variables to 1 in the 180 quadrics above, we 
obtain an affine variety isomorphic to the affine plane. 



Theorem 4.4. The coherent component of the toric Hilbert scheme HilbA 
is isomorphic to ProjC[zG]. 

We present a sketch of the proof of this theorem. The first step is to 
define a morphism from HilbA to P(G) = Proj C [2/1, 2/2, • • • , 2/t]- Consider 
any point / on HilbA- The component / n CGj of I in the i-th Graver 
degree is a linear subspace of codimension 1 in CGj. We can represent this 
hyperplane by a polynomial g%{I) = ^CueG; c « x " which is defined uniquely 
up to scaling. By taking the product of these polynomials, we obtain a 
unique (up to scaling) polynomial which is supported on G = G1G2 • • • G t - 
The assignment I \—* g\{I)g2{I) ■ ■ ■ gt{I) defines a morphism from HilbA 
to P(G). Now consider the restriction of this morphism to the coherent 
component. It is clearly an isomorphism on the (C*) n -orbit of the toric ideal 
I A- To show that it is an isomorphism on the entire coherent component, we 
consider the affine chart around an initial monomial ideal M = in w (lA)- It 
corresponds to the vertex in direction w of the Minkowski sum of the Newton 
polytopes of the polynomials g\{I a) ■, g2{I a) ■, ■ ■ ■ ,gt(lA)- The normal cone 
at that vertex coincides with the cone a of the Grobner fan which has 
w in its interior §3]. The restriction of our morphism to the affine chart 
around M of the coherent component, as described in equation (B), shows 
that this restriction is an isomorphism. The reason for this is that each pair 
of vectors {uj, v{\ seen in the reduced Grobner basis lies in one of the Graver 



fibers Gj. This concludes our sketch of the proof of Theorem 4.4. □ 
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5. Appendix 1: Fourier-Motzkin Elimination 

We now give the Macaulay 2 code for converting the generator /inequality 
representation of a rational convex polyhedron to the other. It is based on 
the Fourier-Motzkin elimination procedure for eliminating a variable from a 
system of inequalities [p3fl . This code was written by Greg Smith. 

Given any cone C C M d , the polar cone of C is defined to be 

C7 V = {x G R d | x • y < 0, for all y G C}. 

For a d x n matrix Z, define cone(Z) = {Zx \ x G ^>q} C K d , and 
affine(Z) = {Zx \ x G M. n } C For two integer matrices Z and -ff, 
both having d rows, polarCone (Z ,H) returns a list of two integer matrices 
{A,E} such that 

cone(Z) + affine(H) = {x G R d \ A l x < Q,E l x = 0}. 

Equivalently, (cone(Z) + affine{H)) y = cone(A) + affine(E). 

We now describe each routine in the package polarCone .m2. We have 
simplified the code for readability, sometimes at the cost of efficiency. We 
start with three simple subroutines: primitive, toZZ, and rotateMatrix. 

The routine primitive takes a list of integers L, and divides each element 
of this list by their greatest common denominator. 

185 : primitive = (L) -> ( 

n := #L-1; g := L#n; 

while n > do (n = n-1; g = gcd(g, L#n) ; 

if g === 1 then n = 0) ; 
if g === 1 then L else apply(L, i -> i // g)); 

The routine toZZ converts a list of rational numbers to a list of integers, 
by multiplying by their common denominator. 

186 : toZZ = (L) -> ( 

d := apply(L, e -> denominator e) ; 
R := ring d#0; 1 := 1_R; 

scan(d, i -> (1 = (l*i // gcd(l,i)))); 
apply(L, e -> (numerator (l*e) ))) ; 

The routine rotateMatrix is a kind of transpose. Its input is a matrix, 
and its output is a matrix of the same shape as the transpose. It places the 
matrix in the form so that in the routine polarCone, computing a Grobner 
basis will do the Gaussian elimination that is needed. 

187 : rotateMatrix := (M) -> ( 

r := rank source M; c := rank target M; 

matrix table(r, c, (i,j) -> M_(c-j-l, r-i-1))); 

The procedure of Fourier-Motzkin elimination as presented by Ziegler in 



[23 1 is used, together with some heuristics which he presents as exercises. 



The following, which is a kind of S-pair criterion for inequalities, comes from 



exercise 2.15(i) in [23]. 



The routine isRedundant determines if a row vector (inequality) is redun- 
dant. Its input argument V is the same input that is used in f ourierMotzkin: 
it is a list of sets of integers. Each entry contains indices of the original rays 
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which do not vanish at the corresponding row vector, vert is a set of in- 
tegers; the original rays for the row vector in question. A boolean value is 
returned. 

i88 : isRedundant = (V, vert) -> ( 

— the row vector is redundant iff 'vert' contains an 

— entry in 'V . 
x := 0; k := 0; 

numRow := #V; — equals the number of inequalities 
while x < 1 and k < numRow do ( 

if isSubset(V#k, vert) then x = x+1; 

k = k+1;) ; 
x 1); 

The main work horse of polarCone .m2 is the subroutine f ourierMotzkin 
which eliminates the first variable in the inequalities A using the double 
description version of Fourier-Motzkin elimination. The set A is a list of 
lists of integers, each entry corresponding to a row vector in the system of 
inequalities. The argument V is a list of sets of integers. Each entry contains 
the indices of the original rays which do not vanish at the corresponding row 
vector in A. Note that this set is the complement of the set Vi appearing in 



exercise 2.15 in [23|. The argument spot is the integer index of the variable 
being eliminated. 

The routine returns a list {proj A,projV} where projA is a list of lists of 
integers. Each entry corresponds to a row vector in the projected system of 
inequalities. The list projV is a list of sets of integers. Each entry contains 
indices of the original rays which do not vanish at the corresponding row 
vector in proj A. 

i89 : f ourierMotzkin := (A, V, spot) -> ( 

— initializing local variables 

numRow := #A; — equal to the length of V 

numCol := #(A#0) ; pos := {}; 

neg := {}; proj A := {}; 

proiV := {}; k := 0; 

— divide the inequalities into three groups, 
while k < numRow do ( 

if A#k#0 < then neg = append (neg, k) 
else if A#k#0 > then pos = append (pos, k) 
else (projA = append(projA, A#k) ; 

proiV = append(projV, V#k) ; ) ; 
k = k+1;); 

— generate new irredundant inequalities. 
scan(pos, i -> scan(neg, j -> (vert := V#i + V#j ; 

if not isRedundant (projV, vert) 
then (iRow := A#i; jRow := A#j ; 

iCoeff := - jRow#0; 

jCoeff := iRow#0; 

a := iCoeff*iRow + jCoeff*jRow; 

projA = append(projA, a); 

projV = append (projV, vert););))); 

— don't forget the implicit inequalities '-t <= 0'. 
scan (pos, i -> (vert := V#i + set{spot}; 

if not isRedundant (proj V, vert) then ( 
projA = append(projA, A#i) ; 
projV = append(projV, vert););)); 

— remove the first column 

projA = apply(projA, e -> e_{l . . (numCol-1) }) ; 
{projA, projV}); 

As mentioned above, polarCone takes two matrices Z, H, both hav- 
ing d rows, and outputs a pair of matrices A, E such that (cone(Z) + 
affine{H)Y = cone(A) + affine{E). 



i90 : polarCone (Matrix, Matrix) := (Z, H) -> ( 
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R := ring source Z; 

if R =!= ring source H then error ("polarCone: " I 

"expected matrices over the same ring") ; 
if rank target Z =!= rank target H then error ( 

"polarCone: expected matrices to have the " I 

"same number of rows") ; 
if (R =!= ZZ) then error ("polarCone: expected " I 

"matrices over 'ZZ'"); 

— expressing 'cone(Y)+aff ine(B) ' as '{x : Ax <= 0}' 

Y := substituted, QQ) ; B := substituted!, QQ) ; 
if rank source B > then Y = Y I B I -B; 

n := rank source Y; d := rank target Y; 

A := Y I -id_(Qq"d); 

— computing the row echelon form of 'A' 
A = gens gb rotateMatrix A; 

L := rotateMatrix leadTerm A; 
A = rotateMatrix A; 

— find pivots 

numRow = rank target A; — numRow <= d 

i := 0; pivotCol := {}; 

while i < numRow do (j := 0; 

while j < n+d and L_(i,j) =!= 1_QQ do j = 

pivotCol = append (pivotCol, j); 

i = i+lj) j 

— computing the row-reduced echelon form of 'A' 
A = ((submatrix(A, pivotCol)) ~ (-1) ) * A; 

— converting 'A' into a list of integer row vectors 
A = entries A; 

A = apply(A, e -> primitive toZZ e) ; 

— creating the vertex list 'V for double description 

— and listing the variables 'T' which remain to be 

— eliminated 

V := {}; T := toList (0 . . (n-1) ) ; 
scan (pivotCol , e -> (if e < n then (T = delete (e, T) ; 

V = append (V, set{e});))); 

— separating inequalities 'A' and equalities 'E' 
eqnRow := {}; ineqnRow := {}; 
scan(numRow, i -> (if pivotCol#i >= n then 

eqnRow = append (eqnRow, i) 

else ineqnRow = append (ineqnRow, i);)); 
E := apply (eqnRow, i -> A#i) ; 
E = apply(E, e -> e_{n. . (n+d-1)}) ; 
A = apply ( ineqnRow , i -> A#i) ; 

A = apply(A, e -> e_(T I toList (n. . (n+d-1) ))) ; 

— successive projections eliminate the variables 'T' . 
if A =!= {} then scan(T, t -> ( 

D := f ourierMotzkin(A, V, t) ; 
A = D#0; V = D#l;)) ; 

— output formating 

A = apply(A, e -> primitive e) ; 

if A === {} then A = map(ZZ"d, ZZ~0, 0) 

else A = transpose matrix A; 

if E === {} then E = map(ZZ"d, ZZ~0, 0) 

else E = transpose matrix E; 

(A, E)); 

If the input matrix H has no columns, it can be omitted. A sequence of 
two matrices is returned, as above. 

191 : polarCone (Matrix) := (Z) -> ( 

polarCone(Z, map(ZZ"(rank target Z) , ZZ~0, 0))); 

As a simple example, consider the permutahedron in M 3 whose vertices 
are the following six points. 

192 : H = transpose matrix{ 

{1,2,3}, 
{1,3,2}, 
{2,1,3}, 
{2,3,1}, 
{3,1,2}, 
{3,2,1}}; 

3 6 
o92 : Matrix ZZ < ZZ 

The inequality representation of the permutahedron is obtained by calling 
polarCone on H: the facet normals of the polytope are the columns of the 
matrix in the first argument of the output. The second argument is trivial 
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since our input is a polytope and hence there are is no non-trivial affine 
space contained in it. If we call polarCone on the output, we will get back 
H as expected. 

193 : P = polarCone (H) 

o93 = (I 1 1 1 -1-1-5 1,0) 
1-11 -5 1 -11 I 
1-1-5 1 -11 1 I 

093 : Sequence 

194 : Q = polarCone (P_0) 

094 = ( I 1 1 2 2 3 3 I , 0) 

12 3 13 12 1 

13 2 3 1 2 11 

o94 : Sequence 

6. Appendix 2: Minimal presentation of rings 

Throughout this chapter, we have used on several occasions the simple, 
yet useful subroutine removeRedundantVariables. In this appendix, we 
present Macaulay 2 code for this routine, which is the main ingredient for 
finding minimal presentations of quotients of polynomial rings. Our code for 
this routine is a somewhat simplified, but less efficient version of a routine 
in the Macaulay 2 package, minPres .m2, written by Amelia Taylor. 

The routine removeRedundantVariables takes as input an ideal I in a 
polynomial ring A. It returns a ring map F from A to itself, which sends re- 
dundant variables to polynomials in the non-redundant variables, and sends 
non-redundant variables to themselves. For example, 

195 : A = QQ[a. . e] ; 

196 : I = ideal (a-b~2-l, b-c"2, c-d~2, a"2-e"2) 

2 2 2 2 2 

o96 = ideal (-b + a - 1, - c + b, - d +c,a -e) 

096 : Ideal of A 

197 : F = removeRedundantVariables I 

8 4 2 

097 = map(A,A,{d + 1, d , d , d, e}) 

097 : RingMap A < A 

The non-redundant variables are d and e. The image of / under F gives 
the elements in this smaller set of variables. We take the ideal of a Grobner 
basis of the image: 

198 : II = ideal gens gb(F I) 

16 8 2 

098 = ideal (d + 2d - e +1) 

098 : Ideal of A 

The original ideal can be written in a cleaner way as 

199 : ideal compress (F. matrix - vars A) + II 

8 4 2 16 8 2 

099 = ideal (d - a + 1, d - b, d - c, d + 2d - e +1) 

o99 : Ideal of A 
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Let us now describe the M&caul&y 2 code. The subroutine f indRedundant 
takes a polynomial /, and finds a variable Xi in the ring of /, such that 
/ = cxi + g, for a non-zero constant c, and a polynomial g which does 
not involve the variable a?j. If there is no such variable, null is returned. 
Otherwise, if Xi is the first such variable , the list {i,c~ l g} is returned. 

1100 : f indRedundant=(f )->( 

A := ring(f ) ; 

p := first entries contract (vars A,f); 

i := position(p, g -> g != and first degree g === 0); 
if i === null then 

null 
else ( 

v := A_i; 

{i;(-l)*(c-(-l)*(f-c*v))} 
) 

); 

The main function removeRedundantVariables requires an ideal in a 
polynomial ring (not a quotient ring) as input. The internal routine f indnext 
finds the first entry of the (one row) matrix M which contains a redundancy. 
This redundancy is used to modify the list xmap, which contains the im- 
ages of the redundant variables. The matrix M, and the list xmap are both 
updated, and then we continue to look for more redundancies. 

1101 : removeRedundantVariables = (I) -> ( 

A : = ring I ; 

xmap := new MutableList from gens A; 
M : = gens I ; 
f indnext := () -> ( 
p := null; 
next := 0; 
done := false; 
ngens := numgens source M; 
while next < ngens and not done do ( 
p = f indRedundant (M_ (0, next) ) ; 
if p =!= null then 

done = true 
else next=next+l; 

); 

p); 

p : = f indnext () ; 
while p =!= null do ( 

xmap#(p#0) = p#l; 

Fl := map(A, A.toList xmap); 

F2 := map(A,A, Fl (Fl .matrix) ) ; 

xmap = new MutableList from first entries F2. matrix; 
M = compress(F2 M) ; 
p = f indnext ( ) ; 

map(A, A,toList xmap)); 
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